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Abstract 

Using both the matrix model prescription and the strong-coupling approach, we describe the 
intersections of n = and n = 1 non-degenerated branches for quartic (polynomial of adjoint 
matter) tree- level superpotential in J\f = 1 supersymmetric S0(N)/USp(2N) gauge theories 
with massless flavors. We also apply the method to the degenerated branch. The general matrix 
model curve on the two cases we obtain is valid for arbitrary iV and extends the previous work 
from strong-coupling approach. For SO(N) gauge theory with equal massive flavors, we also 
obtain the matrix model curve on the degenerated branch for arbitrary N. Finally we discuss 
on the intersections of n = and n = 1 non-degenerated branches for equal massive flavors. 
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1 Introduction 



The exact quantum effective superpotential for the glueball field was proposed by Dijkgraaf 
and Vafa [1, 2, 3] using a zero-dimensional matrix model. Extremization of the effective glueball 
superpotential has led to the quantum vacua of the supersymmetric gauge theory For Af = 1 
supersymmetric U(N) gauge theory with the adjoint matter $, the gauge group U(N) breaks 
into niLi U(Ni) for some n where iV = Yh=i ^i- At low energies, the effective theory becomes 
Af = 1 gauge theory with gauge group U(l) n . The low energy dynamics of this gauge theory 
along the line of [1, 2, 3] have been studied in [4, 5, 6, 7, 8, 9, 10]. 

Although the Af = 2 factorization problem in the strong-coupling approach [11, 12] has been 
solved for small values of iV explicitly [5], the difficulty of the solving the general (arbitrary 
N) factorization problem occurs. In [13], a matrix model curve for U(N) gauge theory with 
cubic tree-level superpotential where the gauge symmetry breaks into U(N) — > U(Ni) x U(N 2 ) 
with N = Ni + N 2 was studied. One of the lessons is that the description of minimization of 
glueball superpotential is useful to obtain a matrix model curve at n — 1 and n = 2 singularity 
2 with arbitrary N for U(N) gauge theory. On the singularity, the two-branch cuts of the 
matrix model curve meet each other. In general, a matrix model curve can be represented 
as = W 2 (x) 2 + fi(x) = Yl i=1 (x — xf) in terms of either a parameter of superpotential and 
fields or the roots of y m . The period integrals on the Riemann surface can be written by the 
elliptic integrals. The glueball equations of motion allowed us to solve for the roots Xi in terms 
of parameters Ni,N 2 and a parameter of superpotential. It was not transparent to obtain the 
four different general roots x,. However at the singularity, since the root x 2 approaches the root 
£3 (reducing of the roots of y m ), by tuning the parameter of superpotential to special values 
depending on Aq and N 2 , it was possible to obtain the roots and parameter of superpotential for 
general (Ni,N 2 ). This overcomes the difficulty in the strong-coupling approach (factorization 
problem) when N is large because a matrix model curve is valid for arbitrary N. 

The extension [5] to the Af = 1 supersymmetric gauge theories with the gauge groups 
SO(N)/USp(2N) was obtained in [14]. The factorization problem gave an explicit construc- 
tions for the matrix model curve for small values of N. In [15], the description for U(N) gauge 
theory was generalized to the SO(N)/USp(2N) gauge theories with quartic tree-level super- 
potential (n — 1) by following the method of [13]. In particular, a new result for degenerated 
branch 3 was obtained besides a simple generalization of the results in [13] to these gauge 

2 For U(N) gauge theory, the tree-level superpotential is written as W($) = Yffk=o ^j-Tr<f> fe+1 . Thus, the 
cubic case can be described by the value n = 2. On the other hand, for SO(N) /U Sp(2N) gauge theories, 
the tree-level superpotential can be written as W(<S>) = Yfltl 

^Tr$ 2A: . The quartic tree-level superpotential, 
which we will deal with mainly, is represented by the value n = 1. 

3 For SO{N)/USp{2N) gauge theories, since a matrix model has Z 2 identification the curve can be written 
in terms of a function of x 2 . One can observe that the position of the branch cut on the fixed point cannot 
move freely by changing a scale A. If there are no branch cuts on the fixed point, called degenerated case, all 



1 



groups. For the degenerated case, the number of roots of y m is less than the number of roots on 
the non-degenerated case. The matrix model curve for the degenerated case 4 can be written 

as 



in terms of either a parameter of superpotential and field or the roots of y m . Using a matrix 



This situation is quite similar to the one studied in [13]. Therefore, the claim in [15] was 
that the glueball approach is suitable for computation of a general matrix model curve on 
the degenerated case which depends on the field F, N and Ni for general breaking pattern 
SO(N) -> SO (N ) x C/(iVi) with N = N + 2N ± . The N can be any generic value. We should 
not restrict our discussion to a special point on a matrix model curve but it covers the whole 
degenerated branch. Thus, as studied in [14] there exist some smooth interpolations between 
the vacua with different breaking patterns like as Coulomb branch in [5]. 

It is natural and interesting to describe the glueball approach by adding the massive or 
massless flavors into the pure SO(N)/USp(2N) gauge theories. Recall that the Seiberg-Witten 
(SW) curves [16, 17, 18] for SO(N)/USp(2N) gauge theories with flavors are characterized by 



where e = for N odd and e = 1 for N even. When the number of flavors Nf becomes zero in 

the power of A and the expressions containing the product for the mass m/ become 1, then the 

above formula reduces to the SW curve of pure gauge theory. Compared with the pure case, the 

factorization problem in the strong-coupling approach with flavors becomes more complicated 

due to the presence of flavor- dependent part. In the strong-coupling approach the Af = 1 

matrix model curve is described by the single root part in the SW curve. One of our aims is 

to observe the singularity on the matrix model curve (i.e., double roots part) for arbitrary N in 

our gauge theories [19, 20] from the point of view of glueball approach. We expect to have a 

the branch cuts move freely on the Riemann surface as illustrated in Fig. 4 in [5]. Therefore, the physics of 
these two cases are different from each other. Thus, we have to study these two cases separately. However, for 
U(N) gauge theory since there is no fixed point, we do not need to study them separately. This is one of the 
reasons why we do not study U(N) gauge theory with flavors in this paper. It would be interesting to study 
U(N) theory flavors in the future. 

4 For the cases of N f = N - 3 for SO(N) and N f = 2N + 1 for USp(2N) gauge theories, the extra piece 
proportional to x 2 should be added. See sections 3 and 5. 




model curve for non-degenerated case, y^ = Yli=o (x 2 + x f )> the degenerated curve can be 
related to y^ = x 1 y 2 md implying that one of the roots in the non-degenerated case vanishes. 
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general matrix model curve for massless flavors without any difficulty because the analysis from 
the strong-coupling approach implies that the SW curve looks like the curve for pure gauge 
theories with different power behavior of x above. However, for massive case (even equal mass 
for flavors), since we have an extra mass parameter in our factorization problem, most of the 
matrix model curves for given N and Nf (the number r in the r-th vacua) in our gauge theories 
do not have an extra double root and therefore, they do not exist at the singularity we are 
interested in. Instead, when they have an extra overall factor x 2 in the matrix model curve 
(i.e., degenerated case), there is a chance to obtain the general matrix model curve because the 
number of roots is reduced and this fact will make easier to compute the equation of motion 
for glueball field. 

In this paper, we deal with two topics: 

• One is n = and n = 1 singularity on the non-degenerated case. This is a simple 
generalization of the results in [13] to the gauge theories with flavors in the vector representation 
for SO(N) gauge theory or fundamental representation for USp(2N) gauge theory. 

• The other is a study for a general matrix model curve for the gauge theories on the whole 
degenerated branch. 

The effective superpotential for the gauge theories with some flavors was already discussed 
in [6, 21, 22, 23, 24]. The contributions of the flavors to the effective superpotential, T 'flavor, 
can be expressed by the integral of matrix model curve (2.1). In [25, 19, 20], the mass of 
flavors possesses the same value and there exists one constraint, W^(±mf) = 0. For the 
SO(N)/USp(2N) cases, since the superpotential W±{x) is an even polynomial in x, the relation 
W^(x) = always has a zero x = as a root. In other words, for massless case, the above 
constraint is automatically satisfied. Therefore, the massless case and the massive case are 
quite different from each other (One cannot obtain the massless case as we take the zero limit 
of the flavor mass in the massive case). As discussed in [19, 20], the massless flavors are charged 
under the factor SO(N ), while the massive flavors are charged under the factor U(Ni) in the 
breaking pattern SO(N) -> SO(N ) x [/(A^). 5 

The organization of this paper is as follows. In section 2, we study the matrix model curve at 
n = and n — 1 singularity for SO(N) gauge theory with massless flavors. This is the simplest 
example, because the contributions coming from the flavors can be written in terms of the dual 
period integral Il . Therefore, using previous results obtained in [15] for the pure gauge theory, 
we obtain a general matrix model curve (2.2) at the singularity for the massless flavor case by 
minimizing the effective superpotential with respect to the glueball field. To demonstrate on the 
validity of general matrix model curve, we deal with some explicit examples for SO(N) gauge 
theories where N = 4,5,6,7,8 or 9 with Nf(< N — 3) flavors. We have checked the precise 

5 Then, we use SO(N ) or U(Ni) to represent the flavors which are charged under the corresponding gauge 
groups. Similar arguments hold for USp(2N) gauge theory. 
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values iV and Ni explicitly by performing the period of T(x) (characterized by a characteristic 
function) over the compact cycles, for some particular examples. The addition map is used. 

In section 3, we study a general matrix model curve on degenerated branch for SO(N) gauge 
theory with massless flavors. Contrary to the results in section 2, we obtain the matrix model 
curve characterized by a single field. There are two classes of matrix model curves, (3.1) and 
(3.2), depending on whether the quantity (N — Nf — 2) is equal to zero or not. After we obtain 
the general matrix model curve for arbitrary N for SO(N) gauge theory, we explicitly check 
the equivalence between two approaches (glueball approach and the factorization problem of 
SW curve). The subtlety for the matrix model curve is discussed when the number of flavor 
Nf reaches the maximum value (N — 3) in the asymptotic region of the theory. We predict 
the matrix model curve for the degenerated SO(N) gauge theories with Nj(< N — 3) massless 
flavors where N >7 and SO(N) -> SO(N ) x [/(A^), N ^ 0. 

In sections 4 and 5, we extend the studies given in sections 2 and 3 for SO(N) gauge theories 
to the USp(2N) gauge theory with massless flavors. The procedures are similar to the SO(N) 
gauge theory. However, there exists one difference: the matrix model curves for the two gauge 
theories can be represented as = x 2 (x 2 ± m) 2 + f 2 (x) where the plus sign corresponds to the 
SO(N) gauge theory and the minus sign corresponds to the USp(2N) gauge theory. Therefore, 
by replacing (N — 2) — > (2iV + 2) and x 2 — > — x\ , we obtain the general matrix model curves 
(4.1), (5.1), and (5.2) for USp{2N) theory from SO(N) theory. When the number of flavors N f 
reaches its maximum value (2iV +1) in a theory, the analysis for the curve in degenerated case 
needs to have an extra term due to the flavor dependent term. We predict the matrix model 
curve for the degenerated USp(2N) gauge theories with Nj(< 2N + 1) massless flavors where 
N > 4 and USp{2N) -> USp(2N ) x U(N 1 ), N ^ 0. 

In section 6, we study the degenerated SO(N) gauge theories with massive flavors. Contrary 
to the massless cases, the constraint for the mass of flavors, W^(±rrif) = 0, is not automatically 
satisfied. However, this condition provides that the contribution arising from flavors can be 
written in terms of a dual period By using this fact and performing the equation of motion 
for glueball field, we obtain a general matrix curve for SO(N) theory. Moreover, we discuss 
r-Higgs branch discussed in [26, 27, 28, 19, 20, 6] in the context of glueball approach. To obtain 
a general matrix model curve for this branch we follow the discussion given in [6] about contour 
of integrals under the change of mass parameter and finally obtain the effective superpotential 
for the branch. Minimizing the effective superpotential we derive the matrix model curve for the 
theories (6.5), (6.6), (6.8), and (6.9). After the general discussions, we deal with some explicit 
examples, SO(N) with N = 4, 5, 6 and check the explicit agreement between two approaches. 
We predict the solutions for the degenerated SO(N) gauge theories with Nf(< N — 3) massive 
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flavors where N > 7 and SO(N) -> S'O(JVo) x ^(iVi), N + 0. 6 

In section 7, we give some comments on the n = and n — 1 singularities for massive 
SO(N)/USp(2N) cases. We do not obtain the general matrix model curves because, in general, 
the J 7 flavor cannot be written as a dual period ITi. Therefore, the matrix model curve cannot be 
represented by any simple formulation. However, when the condition, c < where c is defined 
in (7.10) or (7.11), should be satisfied, the flavor-dependent part can be written in terms of ITi. 
Thus, in this particular case, we obtain a matrix model curve for arbitrary N. 

There exist many related works from the different matter representations for various gauge 
theories [29]- [37]. 

2 The n = and n = 1 singularity: SO(N) gauge theory 
with massless flavors 

Let us consider SO(N) gauge theory with quartic tree- level superpotential (n — 1) and Nf 
flavors. We are interested in the intersections of the n = and n — 1 branches. At these special 
points, a vacuum with unbroken gauge group SO(N ) x U(Ni) meets a vacuum with unbroken 
gauge group SO(N) with N = N + 2N X . Since we are considering an asymptotic free theory, 
we restrict the number of flavors Nf to satisfy the condition Nf < N — 2. As already clarified 
in [19], the relation between W%{x) = x (x 2 + m) and F 6 (x) corresponding to an M = 1 matrix 
model curve is characterized by 

y 2 m = F 6 (x) = W^x) 2 + f 2 (x), N,<N-3, 
y 2 m = F 6 (x) = W^(x) 2 + f 2 (x)-Ax 4 A 2 , N f = N - 3 

where /^(a?) = f2^ 2 + /o- These relations are valid for any SO(N) gauge theories where N is 
even or odd. All masses of the flavors vanish. The A is an M = 2 strong-coupling scale. The 
fluctuating fields f 2 and fo relate to the glueball fields Si and So respectively. 

The effective superpotential for the SO(N) gauge theory with Nf flavors is given by [22, 
21, 23, 24] 

W cS = 2m [{N - 2)II + 2^11!] - 2(N - 2 - N f )S log (A) _ j=f lavor , 
where the contributions from the flavors are given by the integral for y m over x: 

Mo 

^flavor = 2^ Vmdx. (2.1) 

f=l J±m f 

6 As one can see in [20], there is no degenerated branch for USp(2N) theory. Therefore, we do not need to 
study them here. 
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The log-divergent term was renormalized by the bare Yang-Mills coupling and the cut-off A was 
replaced by the physical scale A. 7 Note that there exists a relation S — S + 2Si. The periods 
Si and dual periods IT of holomorphic 3-form for the deformed geometry were written as the 
integrals over the x-plane [15, 38, 39]. The m/ is the mass of the f-th. flavor and we have chosen 
the location of holomorphic 2-cycles at x — ±m,f. Around these points, D5-branes are wrapping 
on holomorphic 2-cycles, providing a source for a single unit of RR flux at x — ±m/. Thus the 
integral of 3-form RR flux around x = ±rrif is equal to 1. The full effective superpotential can 
be generated by adding this contribution to the effective superpotential for pure gauge theory. 

First, let us consider the simple case in which all masses of flavors are zero in this section. 
Later, in sections 6 and 7, we will deal with the massive cases. Since the contribution from the 
flavors, Flavor, becomes 2itiNfT\.Q ) we can use the results for pure SO(N) gauge theory studied 
in [15] without much efforts. We only have to replace the quantity (JV — 2) appeared in the 
matrix model curve in [15] with (N — Nf — 2). Therefore, by combining the flavor-dependent 
part with the dual period il term, the effective glueball superpotential for massless flavors is 

W cS = 2m[(N -2-N f )U + 2N 1 U 1 ]-2(N-2-N f )S\og^y 

The periods of T(x) defined as the periods of the 1-form at the n = and n — 1 singularity are 
given similarly and together with different power behavior of x and A, due to the presence of the 
flavors, the T(x) also takes the form as a complete differential. With the above replacement 
in mind, we immediately obtain the matrix model curve near the singularity on the non- 
degenerated case from the glueball equations of motion: 



yl = x 2 (x 2 + m) + (f 2 )x 2 + </ > = [x 2 + 2 V A 2 (1 + c)] (x 2 + A V A 2 ) , 

where the fields, parameter of a superpotential, and the glueball field are given by 

(/ 2 ) = VA 4 (l + 2c), (/o) = 16^ 3 A 6 (l + c) 2 , 
m = 2 V A 2 (2 + c), (S) = - V 2 A A (l + 2c), c = cos ( - - ) . (2.2) 

Here r] is the (N - N f - 2)-th root of unity for (N - N f - 2) even and the (N - N f - 2)-th 
root of minus unity for (N — Nf — 2) odd. That is, 

v N-N f -2 = x for ( iV _ iV/ _2) everi; 

= _ x for ( iV _ iV/ _2) dd. 



7 For the tree-level superpotential of degree 2(n + 1), the effective superpotential for the SO(N) gauge theory 
with N f flavors is W cB = 2m [{N - 2)U + J27=i 2N i U i\ - 2 (JV - 2 - JV» (S + E^i lo S (^) - ^flavor- 
For quartic case (n = 1), this will lead to the above effective superpotential. 
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This matrix model curve with massless flavors is exactly the same as pure case with the replace- 
ment (N — 2) — > (N — Nf — 2). 8 Note that the above solutions have the special combination 
(N — Nj — 2) which reflects the addition map discussed in [19], relating the Chebyshev branches 
(or the Special branches) of two different gauge groups. Namely, we have seen that the Cheby- 
shev branch of SO(N) with Nf massless flavors can be reduced to the Chebyshev branch of 
SO(N') with Nf massless flavors when there exists a relation: N — Nf — 2 = N' — Nf — 2. 

Let us demonstrate these general features (2.2) explicitly by comparing them with the results 
from strong-coupling approach obtained in [19] already. Let us consider SO(N) gauge theories 
with JV/(< N-S) flavors where TV = 4, 5, 6, 7, 8 or 9. 

• SO (A) with N f = 1 

In [19], the factorization problem resulted in the matrix model curve y^ 2 which contains an 
overall factor x 2 . 9 Therefore, there exists no solution for non-degenerated case. We will come 
to this for degenerated case later. 

• 50(5) with N f = 1 

In [19], the matrix model curve y^ 2 contains an overall factor x 2 . Therefore, there exists no 
solution for non-degenerated case. We will come to this for degenerated case later. 

• 50 (5) with N f = 2 

The SW curve is the same as the curve of SO [A] with Nf — 1 which reflects the addition 
map. Then the factorization problem will lead to the one in 5*0(4) with Nf = 1 theory. 

• 50(6) with N f = 1 

This example is the first nontrivial case. By putting the values (N, N±, Nf) = (6, 1, 1) into 
the curve (2.2), we predict the matrix model curve for this case, 

y 2 m = x 2 [x 2 + 3r?A 2 ) 2 - 4A 6 = (x 2 + r]A 2 ) 2 (x 2 + 4r?A 2 ) 

where r\ satisfies rf = — 1 or equivalently rf = 1. Recall that the matrix model curve for pure 
5*0(5) gauge theory [15] is identical to this curve. That is, N — Nf = 5. The factorization 
problem [19] resulted in the matrix model curve y^ 2 = x 2 (x 2 — a 2 ) - 4A 6 . There exists 
a symmetry breaking 50(6) — > 50(4) x U(l). For a 2 = —3r]A 2 , this can be factorized as 
(x 2 + Ar]A 2 ) (x 2 + r]A 2 f. For a 2 = 2A 2 , the curve x 2 (x 2 - a 2 ) 2 - a 4 x 2 will be x 4 (x 2 - 4A 2 ). 
The latter implies the degenerated case. The characteristic function Pq{x) = x 4 (x 2 + 3??A 2 ) 
becomes 2x 3 p 3 A 3 % (^a)- These are the vacua surviving when the N = 2 theory is perturbed 
by a quadratic (n = 0) superpotential and the 50(6) gauge theory becomes massive at low 
energies. 

8 0nc expects that the matrix model curve for U(N) gauge theory with massless flavors (U(N) — > 
U(Ni) x £/(iV 2 ) where N = N\ + N 2 ) at n = 1 and n = 2 singularity takes the form = 

[x + r]A(2 + c)] [a; - r?A(2 - c)] {x - r]Ac) 2 where c = cos (wBij) and ^ = L 

9 We denote the curve from strong-coupling approach (factorization problem) by y^ 2 and the curve from 
glucball approach by y^- 
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Let us demonstrate how one can obtain the precise values for (N ,Ni). The results for a 
matrix model curve and a characteristic function can be rewritten as follows (we put A = 1 
and 77 = — 1) : 



&n 2 = (x 2 -A)(x 2 -l) , P 6 (x) = X A (x 2 - 3) . 

By using these relations, the function T(x) = Tr^-^ is given as [19] 

T(x) - n{x) ~ iJ ^ P *W 1 ^/ + 2 ) _ 

^Pe(x) 2 -4x 6 A 6 x ^x 2 -4 x 

There exist three branch cuts on the x-plane [—2, —1 — e], [—1, 1], and [1 + e, 2] before taking 
the limit x\ — > x (e — > 0). 10 Since we are assuming n = and n = 1 singular case, these 
three-branch cuts are joined at the locations of x — ±1 after taking the limit and they reduce 
to a single-branch cut [—2,2]. Therefore, we can explicitly calculate the values (N ,Ni) as 
follows: 

N = — i T(x)dx = — - _ + -)dx=[- n +3 = 4, 



2iri Ja 2ni J- 

2 

2tH 



Ni = — I T(x)dx 




as we expected. 11 

• SO(Q) with Nf = 2 

Once again the factorization problem [19] resulted in the matrix model curve that contains 
an overall x 2 factor implying a degenerated case. 

• SO (6) with N f = 3 

The factorization problem [19] turned out the matrix model curve that has an overall x 2 
factor and is classified as a degenerated case which will be studied later. 

• SO(7) with N f > 1 

The curve for SO(7) with Nf flavors is the same as 5*0(6) curve with (Nf — 1) flavors. In 
particular, since the SW curve of SO (7) theory with Nf = 1 is identical to the 50(6) gauge 
theory with Nf = 0, the factorization problem in [14] implies that the matrix model curve is 

y 2 m = x 2 (x 2 + Ari A 2 ) 2 + Ar] 2 A A x 2 + 16?? 3 A 6 = (x 2 + 2??A 2 ) 2 (x 2 + 4??A 2 ) , = 1 



10 Recall that a matrix model curve is given by — YYi=o { x2 + X X)- 

"Recently [40], the construction of T(x) on an elliptic curve in terms of three discrete parameters and two 
continuous ones was studied in the breaking pattern U(N) — > U (JVi) x U(N-2) where N = N\ + N 2 - This method 
will give some hints to obtain this 1-form explicitly without resort to a characteristic function we did here. 
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where the breaking pattern 5*0(7) — > 50(5) x U(l) exists. 12 At the intersections with n = 
branch, the characteristic function P%(x) can be written in terms of the Chebyshev polynomial 
of degree 4 with an appropriate identification between these two functions. 

• 50(8) with N f > 1 

The curve for 50(8) with Nf = 1 can be obtained from the factorization problem. The result 
is, by applying the method of [19] to the present case here, y^ 2 = {x 2 + ^^A 2 ) (x 2 + 4r]A 2 ) 
with i] 5 = — 1. This curve is exactly the same as the expression from the glueball approach 
by inserting the values (JV , N u N f ) = (6, 1, 1) into (2.2) where 50(8) -> 50(6) x U(l). Note 
that the above matrix model curve is the same as the curve for pure 50(7) gauge theory 
[15]: N — Nf = 7. Although there is no information about the curve with larger flavors from 
factorization problem, one predicts a matrix model curve through (2.2). 

• 50(9) with Nf > 1 

The curve for 50(9) with Nf flavors is the same as 50(8) curve with (Nf — 1) flavors. In 
particular, since the SW curve of 50(9) theory with Nf = 1 is identical to the 50(8) gauge 
theory with Nf = 0, the factorization problem in [14] implies that the matrix model curve 
becomes 

y 2 m = x 2 [x 2 + ^N 2 ) 2 + 8r] 2 A A x 2 + 36r? 3 A 6 = (x 2 + 3??A 2 ) 2 (x 2 + 4r?A 2 ) , if = 1 

where the breaking pattern 50(9) — > 50(7) x U(l) exists. 13 At the intersections with n = 
branch, the characteristic function Ps{x) can be written in terms of the Chebyshev polynomial 
of degree 6 with an appropriate identification between these two functions. These are the vacua 
that survive when the Af = 2 theory is perturbed by a quadratic superpotential (n — 0). 

Recalling the replacement of (N — 2) in pure gauge theory result with (N — Nf — 2) 
for massless flavors we can get the coupling constant near the singularity without explicit 



12 Let us check how one can determine the precise values for (No, N±). The results for a matrix model curve 
and a characteristic function are given by (wc put A = 1 and r\ = —1): y^ 2 = (x 2 — 4) (x 2 — 2) ,P%(x) = 

x 2 (x 4 - 4x 2 + 2). By using these relations, the function T(x) is given by T(x) = P ^" ( ^^ =^ + (JV/ + 2) = 

— ^ x 2_ 4 + f • Note that there was a typo in T(x) [14] and there should be present an extra 1/x term in T(x). 

There exist three-branch cuts on the x-plane [—2, — v/2 — e], [— v/2, y/2], and [\/2 + e, 2] before taking the limit 
Xi — > xo(e — > 0). These three-branch cuts are joined at the locations of x = ±\/2 after taking the limit and 
they reduce to a single-branch cut [—2, 2]. Therefore, we can explicitly compute the values (No, N\) as follows: 

13 By using y^ 2 = (x 2 - 4) (x 2 - 3) 2 ,P s (x) = x* (x 2 - 3) 2 - 2x 2 , the function T(x) is given by T(x) = 
Ps(j) — p p^l + (Ei±2l = _ Thcrc cxigt thrce _ branch cuts on the x _ p i ane [_2, — y/3-e], [-VH, VS], 

and [>/3+e, 2] before taking the limit x\ — > x (e — » 0). These three-branch cuts are joined at the locations of x = 
±V3 after taking the limit and they reduce to a single-branch cut [—2, 2]. Therefore, we can explicitly compute 

the values (N , N,) as follows: N = & (--^_ + §) dx = 7, N, = J% (-^|= + f ) dx = 1. 
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calculation, 

J_ d 2 F p {S k ) 1 " d 2 F p {S k ) ( N Q -N f -2 \ d 2 F p (S k ) 
2m ij dS t dS, j Nif^ 1 OS.dSi ij \ N t J dSodSi ' 

where i, j = 1, 2, • • • , n. In the quartic tree-level superpotential case (n = 1) we are considering, 
there is only one coupling constant and it is given by, due to the cancellation of first two terms 
above, 



1 in(N -N f -2 



,2 



2ni ' 16 (N - N f - 2f log (^) ' 
where k' = — §V-| — \{. Since we are assuming asymptotically free gauge theory, the quantity 

x l \ x 2~ x o) 

(N — Nf — 2) is greater than zero. When the condition (N — Nf — 2) = is satisfied, r becomes 
zero. In this case, the asymptotic freedom for SO(N ) gauge theory breaks down. Therefore 
the situation is the same as the one-cut case (equivalently n = case) in which the coupling 
constant is trivially zero. Our result gives a consistency. 

3 The curve for degenerated SO(N) gauge theory with 
massless flavors 

In the non-degenerated case, every root of W^(x) possess at least one D5 brane wrapping 
around it. For the degenerated case, some roots do not contain wrapping D5 branes around 
them. At the classical limits, there are symmetry breakings characterized by 

SO(N) -> SO(N ) x [/(TVi), SO(N) -> U([N/2}), SO(N) -> SO(N). 

The last pattern occurs only if further constraints are imposed on the first two breaking pat- 
terns. For degenerated case we can write the matrix model curve which is valid on the whole 
degenerated branch. As we did for pure gauge theory, the dual periods are defined by the 
integrals over x and the equation of motion of a single field F provides one relation. Recalling 
the simple replacements (N — 2) -> (N - N f - 2) and (JV - 2) -> (N - N f - 2) we made 
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before, we can obtain the matrix model curve as 



2 / 2 \ 2 X 2 -F 



:-f) 



N -Nf-2 



(3.1) 



Here we rescaled the function and corrected a typo appeared in [15]. This formula depends 
on the parameter F, N , N 1: and Nf. Turning on a field F to the particular values will reduce 
to the last breaking pattern SO(N) -> 5O(A0 above. When (jV -N f -2) = (we will meet 
this situation in the examples below), the above formula becomes ill-defined. We go back to 
the derivation of the formula 15 and can find another formula for this special case, 

y 2 m , d =(x 2 + D) 2 -^ ^ = 1, for (JV -JV)-2) = (3.2) 

where D is arbitrary parameter. 

To demonstrate these general solutions (3.1) and (3.2), let us consider SO(N) gauge theories 
with Nf(< N — 3) flavors where N = 4, 5 or 6. 

• SO (A) with Nf = 1 

There are two breaking patterns on the degenerated case: 5*0(2) x U(l) and U(2). 

1) For the breaking pattern 50(4) -> 50(2) x U(l), by substituting the values (N, N u N f ) = 
(4, 1, 1) into (3.1) with the footnote 14 (recall that N = Nq + 2Ni), the matrix model curve 
from glueball approach can be represented by 

y 2 m ,d= (x 2 + m - 2A 2 f + AF, m = J + A 2 . 

We find that the solutions given in [19] perfectly agree with the above result for y 2 nd from 
glueball approach. In other words, the factorization problem characterized by P 2 (x) — 4x 6 A 2 = 
H 2 (x) {x 2 — D) 2 + AF reproduces the above solution where D = A 2 — ^ (For the notation 
of [19], a is equal to -2D and (a 2 + 2A 2 ) corresponds to A 2 - = D). Then a 2 = -A 2 - 
which is equal to —to. For F — 0, there exists a Special branch with 50(4) — > 50(4) and for 
(to - 2A 2 ) 2 + AF = 0, there exists a Chebyshev branch with 50(4) -> 50(4). 
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For the degenerated case also, the relation between W^x) and F^x) can be written as 
vIa = ^ 4 (*)=(^) 2 +4F, N f <N-3, 
V 2 m , d = ^ 4 (x)=('^MV+4F-4, 2 A 2 , 7V / = 7v-3. 



In particular, when Nf = N — 3, the matrix model curve y^ d has an extra contribution from 4x 2 A 2 . Therefore 
in this particular case, the curve should be d = (x 2 + m — 2A 2 ) 2 + 4F and m — 2 A 2 = K ^ F . 

15 In this case, there is a relation 2A^ log | grr^ |=0- By putting 77 = ± fc ^ a ^ and realizing D = | (a 2 + fe 2 ) , 
one gets the above formula. 
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2) For the breaking pattern SO (A) — > U(2), we obtain the following matrix model curve 
from the glueball approach, by plugging the values (N,Ni,Nf) = (4,2, 1) into (3.1) with the 
footnote 14, 



Again by using the different parametrization (as in previous case) from the one in [19], we 
can obtain the same matrix model curve with D which is equal to — m + 2A 2 above where the 
relation 6a 2 A 2 + 14A 4 - 2 (a 2 + 3A 2 ) V4a 2 A 2 + 9A 4 = 47 (corresponding to b - a 2 /4 in the 
notation of [19]) will provide one solution for a 2 = —m in the above. Therefore, the results 
obtained from two approaches are equivalent to each other. 
• 50(5) with Nf = 1 

There are two breaking patterns on the degenerated case: 5*0(3) x U(l) and 5*0(1) x £7(2). 
1) For the former case 50(5) -> 50(3) x U(l), since the condition (JV - N f - 2) = is 
satisfied, the matrix model curve can be read off, from (3.2), 



This is exactly the first kind of solution given in [19] and the D here corresponds to a/2 = —a 2 
there and — 4A 4 corresponds to (b — a 2 /A). For D 2 = 4A 4 , there is a Chebyshev vacuum where 



2) For the second breaking pattern 50(5) — > 50(1) x U(2), by paying attention to the 
value of N , we simply put N — 1 and find the following matrix model curve, 



where r\ is 2-nd root of unity. These solutions are exactly the same as the ones from factorization 
problem in [19] (m = a/2 and b — a 2 /A = 47). For a = ±2A 2 (+ sign corresponds to i] = — 1 
above and — sign corresponds to rj — 1), there is a Special vacuum where 50(5) — > 50(5). 
This curve is the same as the matrix model curve for pure degenerated 50(4) gauge theory 



[15] (N — Nf — 4) where 50(4) -> U(2). 
• 50(6) with N f = 1 

In this case, there are three kinds of breaking patterns on the degenerated case: 50(4) x 
17(1), 50(2) x 17(2), and 17(3). 

1) For the breaking pattern 50(6) — > 50(4) x U(l), the matrix model curve is represented 
as from (3.1) with (77, N u N f ) = (6, 1, 1) 





50(5) -> 50(5). 
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Identifying AF here with b in [19] (and m = —a), we can see the exact agreement. Recall that 
this curve is identical to the matrix model curve for pure degenerated 5*0(5) gauge theory [15] 
(N — Nf — 5) where 50(5) -> 50(3) x U(l). 

2) On the other hand, for the breaking pattern 5*0(6) — > U(3), the matrix model curve can 
be written as from (3.1) together with (N, N 1: Nf) = (6, 3, 1) 

y 2 m ,d= (x 2 + m) 2 + AF, m = HZ-r l 2 A 2 

where if = 1. By identifying (a, b) with (— m, AF), we can see the complete agreement with 
the first, second and third solutions of the equation (6.38) in [19] although we have not written 
them explicitly. For b = 4F = 0, this leads to a Special vacuum of 50(6) — > 50(6). 

3) For the breaking pattern 50(6) — > 50(2) x U(2), the solution from factorization was 
not new one but can be written as the previous solution we have discussed already. 

• 50(6) with Nf = 2 

1) For the breaking pattern 50(6) — > 50(4) x U(l), the matrix model curve can be written 
as from (3.2) 

This corresponds to the first kind solution in [19]. When D 2 = 4A 4 , this will lead to a Chebyshev 
vacuum with 50(6) — > 50(6). 

2) For the breaking pattern 50(6) — > 50(2) x U(2), the matrix model curve is given by 
from (3.1) with (N,N u N f ) = (6,2,2) 

Vm,d= (x 2 + mf + AF, m =HZ- v A 2 

where i] is 2-nd root of unity. This can be seen from the degenerated pure 50(4) gauge theory 
with the breaking pattern 50(4) — > U(2) [15] through the addition map. If we change the 
parametrization by 4F = 4i]aA 2 + 4A 4 , the above matrix model curve can be rewritten as 

y 2 m ,d= (x 2 + af + 4 V A 2 x 2 , 

which is exactly the same as the second kind of solution in [19]. For a = 0, it gives a Chebyshev 
vacuum of 50(6) — > 50(6). For a = —r]A 2 , a Special branch appears. 16 

• 50(6) with N f = 3 

In this case, there are three breaking patterns: 50(4) x U(l), 50(2) x U(2) and U(3). 
16 For the breaking pattern SO(6) — > U(3), we can predict the following matrix model curve, from (3.1) 

3 

with (TV, Ni,Nf) = (6, 3, 2), y 2 m d = (x 2 + m) 2 + AF, and m = ^ - |AFi where rf = 1. In this case, the 
factorization problem becomes very complicated. Since we cannot obtain the explicit results for the matrix 
model curve from the factorization, we have not checked them in this example. 
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1) For the first two cases, we expect that the matrix model curves are the same as the 
ones in SO (A) theory with Nf = 1 as discussed above, which comes from the addition map for 
massless case. For the breaking pattern 5*0(6) — > 50(4) x U(l), the matrix model curve can 
be represented as from (3.1) with the footnote 14 

y 2 m, d = (x 2 + m- 2A 2 ) 2 + AF, m = + A 2 . 

From the result of [19], the matrix model curve was y^ 2 = (x 2 — a 2 — 2A 2 ) 2 — 4a 2 — 4A 4 . By 
introducing a new field AF = —4a 2 A 2 — 4A 4 , one can identify a 2 = —A 2 — which is equal to 
— m in the above. 

2) For the breaking pattern, 50(6) -> 50(2) x U(2), from (3.1) with the footnote 14 the 
matrix model curve can be obtained, with (N, N±, Nf) = (6, 2, 3), 




y 2 m4 =(x 2 + m-2A 2 ) + AF, m = (-—) - F - — + 2A 



As expected, this matrix model curve is the same as the one in 50(4) theory with Nf = 1 
case and agrees with the one from the factorization problem. In other words, the factorization 
problem Pq(x) — Ax 4 A 4 = x 4 H 2 (x) (x 2 — D) 2 + AF reproduces the above solution where D = 

- + F (-fi)*. When F = 0, it gives a Special vacuum of 50(6) -> 50(6). For 

(m - 2A 2 ) 2 + AF = 0, it leads to a Chebyshev vacuum of 50(6) -> 50(6). 17 

In summary, for the SO(N) gauge theories with Nj(< N — 3) flavors where = 4, 5 or 
6, we have convinced ourselves that the solutions (3.1) and (3.2) coincide with the matrix 
model curve from the strong-coupling approach except that we could not check the 50(6) 
gauge theory with Nf = 2, 3 where there exists a breaking pattern 50(6) — > U(3) (in this 
case, A"o = 0). The main obstacle why the factorization problem for SO(N) — > U([N/2\) 
becomes so complicated was the SW curve does not possess any power of x 2 , contrary to the 
factorization problem for 50(A^) -> 5O(iV ) x ^(A^A^o ^ in which the SW curve has a 
power of x 2 (either Chebyshev or Special branch) making the factorization easier. Therefore, 
at least one expects that our matrix model curve (3.1) and (3.2) will predict the solutions for 
the degenerated SO(N) gauge theories with Nf(< N — 3) massless flavors where N > 7 and 
50(AO -> 5O(iV ) x U(N 1 ), N ^ 0. 

17 The last one is the 5*0(6) — ► U(3) breaking pattern. In this case the matrix model curve is given by, 
from (3.1) and the footnote 14 with (N, iVi, N f ) = (6, 3, 3), y 2 m d = (x 2 + m - 2A 2 ) 2 + AF, and m = 5 - 

F (^j^j 5 + 2A 2 . Since the factorization problem becomes very complicated, we have not checked them in this 
example. 
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4 The n = and n = 1 singularity: [/ Sp(2N) gauge theory 
with massless flavors 



As in the SO(N) case discussed previous section, we can obtain the matrix model curve 
without explicit calculation by looking at the pure U Sp(2N) gauge theory. We only have to 
replace the expression (2N + 2) in the result of [15] for the matrix model curve with (2N — 
Nf + 2). With this replacement in mind, we immediately obtain the matrix model curve near 
the singularity on the non-degenerated case from the glueball approach 18 : 

y 2 m = x 2 (x 2 + mf + (f 2 )x 2 + (/o) = [x 2 + 2 V A 2 (1 + c)f (x 2 + 4 V A 2 ) , 

where the fields, parameter of a superpotential, and glueball field are given by 

(/ 2 ) = VA 4 (l + 2c), </„> = 16^ 3 A 6 (1 + cf , 
m = 2 V A 2 (2 + c), (S) = VA 4 (l + 2c), c ^ cos ( _ ^ _ ^ , ( 4 .1) 

Here 77 is the (2iV - N f + 2)-th root of unity, that is rf N - N f +2 = 1. 19 This matrix model 
curve with massless flavors is exactly the same as pure case with the replacement (2N + 2) — > 
(2N — Nf + 2). Note that the above solutions have the special combination (2N — Nf + 2) 
which reflects the addition map discussed in [20], relating the Chebyshev branches (or the 
Special branches) of two different gauge groups. Namely, the curve of USp(2N) with Nf 
massless flavors can be reduced to the curve of U Sp(2N — 2r) with (Nf — 2r) massless flavors. 

Let us demonstrate these general features (4.1) explicitly by comparing them with the results 
from strong-coupling approach obtained in [20] already. Let us consider U Sp(2N) gauge theories 
with jV/(< 2iV + 1) flavors where N = 2 or 3. 

• USp(A) with N f = 1 

From the results in [20], the matrix model curve has an extra double root at d = 

on the non-degenerated case. It was factorized as y^ 2 — (t + ^^eA 2 ) (t + 4eA 2 ) where e is 

5-th root of unity. 20 On the other hand, the matrix model curve from the glueball approach 

18 The effective glueball superpotential for massless flavors is W eS = 2ni [(2N + 2 - N f )U a + 27Vini]-2(2A^+ 
2-AT / )Slog(A). 

19 By recognizing the matrix model curve for USp(2N) gauge theory characterized by = Yli=o { x2 ~ 
one can read off the matrix model curve for USp(2N) with flavors from the curve for SO(N) gauge theory 
with flavors by replacing (N — 2) — > (2iV + 2),xf — > —x\ where i = 0,1,2 as follows: x\ = — 47yA 2 , and 
x\ = -2r]A 2 (1 + c) where 77 satisfies the condition ry 2Ar - w /+ 2 = 1 for (27V - N f + 2) even and ry 2Ar - w /+ 2 = -1 
for (27V — Nf + 2) odd. When we define e = —r\ it satisfies e 2N ~ N f+ 2 — 1 for both cases. The solutions can 
be rewritten as x\ = 4eA 2 , x\ = 2eA 2 (1 + c), and e 2N ~ N f+ 2 = 1 where c is given in (4.1). This is exactly the 
same as (4.1) by using the different parametrization. 

20 We use a relation t = x 2 all the time in this paper. 
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becomes, from (4.1) with (N,N u N f ) = (2, 1, 1), 



2 

where the other values in different parametrization are given by 
(7±V5) 



m 



- V A 2 , (f 2 ) = 2 (l ± y/5) rfk\ (/„) = 2 (7 ± 3Vb) ^ 3 A 6 . 



2 

Here 77 is 5-th root of unity. Now it is easy to see the exact agreement between two results. 
• USp(A) with N f = 2 

The matrix model curve obtained from the factorization problem can be represented, at 
singular point, as y^ 2 — (t + 2i]A 2 ) 2 (t + Ar/A 2 ) where 77 is 4-th root of unity (Note that there 
is a typo in [20]). The characteristic function B e (x) = x 2 (x 4 + Ar\A 2 x 2 + 2?i 2 A 4 ) will become 
2p 2 x 2 A 4 T 2 (^2 + l) with p 4 = 1 by identifying p = rp These are the vacua surviving when the 
Af = 2 theory is perturbed by a quadratic (n = 0) superpotential and the U Sp(A) gauge theory 
becomes massive at low energies. To see the correctness of our formula (4.1) we put simply the 
values (2N , N u N f ) = (2, 1, 2) (Note that 2N = 2N + 2A^) where USp{A) -> USp{2) x 17(1) 
into the formula (4.1) and get the following matrix model curve, 

y 2 m = x 2 (x 2 + 4?pA 2 ) 2 + Ar] 2 A 4 x 2 + 16?? 3 A 6 = (x 2 + 2rp\ 2 ) 2 (x 2 + 4?pA 2 ) . 

Now there exists a perfect agreement. 

Let us demonstrate how one can determine the precise values for (2N 0) Ni). The results can 
be rewritten as follows (we put A = 1 and 77 = 1) : 

y^ 2 = (V + 4) (x 2 + 2) 2 , B 6 (x) = x 2 (x A + 4x 2 + 2 

By using these relations, the function T(x) = Tr^-^ is given as [20] 

T[x) B' (x) - ^B 6 (x) | (N f - 2) 4 



sjB 6 (x) 2 - Ax 4 A 8 % ^/x 2 + A 

There exist three branch cuts on the x-plane [— 2i, —(y/2 + e)i], [— \/2i, \/2i], and [(\/2 + e)i, 2i] 
before taking the limit x\ — > xq. Since we are assuming n = and n = 1 singular case, these 
three-branch cuts are joined at the locations of x — ±y/2i after taking the limit and they reduce 
to a single-branch cut [— 2i,2i\. Therefore, we can explicitly calculate the values (2N ,Ni) as 
follows: 

If 2 A 

2N = <b T(x)dx = / dx = 2, 

2m Jao v 1 2m J-V2i ^x 2 + A 

lr 2 r 2i A 

Ni = <p T(x)dx = / dx = 1 

2m J Al y ' 2m ^/x 2 + A 
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as we expected. 

• USp(A) with N f = 3 

The matrix model curve at the singularity was given in [20] as y m = (t + r]A 2 ) (t + 4t]A 2 ) 
where rf = 1. Note that a 2 corresponds to -3i]A 2 . The B 6 (x) is equal to 2p 3 x 3 A 3 T 3 (^) 
by putting p — rj. At these intersection points, the vacua survive when the J\f = 2 theory is 
perturbed by a quadratic (n = 0) superpotential. If we put the values (2iV , N 1: Nf) = (2, 1, 3) 
where USp(A) — > USp(2) x U(l) into our formula (4.1) we obtain the matrix model curve at 
the singularity, 

y 2 m = x 2 (x 2 + 3??A 2 ) 2 + 4?? 3 A 6 = (x 2 + ??A 2 ) 2 (x 2 + 4??A 2 ) . 

We see the agreement between two approaches. 

• USp(A) with N f = 4 

There is no non-degenerated case and we will discuss it in the degenerated case later. 

• USp(A) with N f = 5 

In [20], the factorization problem resulted in the matrix model curve y^ 2 which contains an 
overall factor x 2 . Therefore, there exists no solution for non-degenerated case. We will come 
to this case for degenerated case later. 

• USp(<6) with Nf = 2 

1) For the first solution USp(Q) — > USp{A) x U(l) in [20], the matrix model curve has 
an extra double root at d = — eA 2 , where e is 6-th root of unity. At the singular points, 
the matrix model can be factorized as y^ 2 = (t — 3eA 2 ) 2 (t — 4eA 2 ). If we simply put the 
values (2N , N±, Nf) = (4, 1, 2) into our general formula (4.1) for the matrix model curve at the 
singularities, we obtain the following results, 

yl = x 2 (x 2 + 5 V A 2 ) 2 + 8 V 2 A 4 x 2 + 36?? 3 A 6 = (x 2 + 3^A 2 ) 2 (x 2 + 4 V A 2 ) 

where r] is 6-th root of unity. If we identify this i] with — e in the solution for factorization 
problem, we see the agreement of two approaches. 21 This choice for B 8 (x) can be written as a 
Chebyshev polynomial of degree 3 with appropriate identification. At these intersection points, 
the vacua survive when the M = 2 theory is perturbed by a quadratic (n = 0) superpotential. 
By the addition map, the above matrix model curve from glueball approach can be obtained 

21 By using y^ 2 = (x 2 + 4) (x 2 + 3) 2 , and B$(x) = x 2 (x 6 + 6x 4 + 9x 2 + 2) , the function T(x) is given 
by T(x) = B 8(*)—frB s (x) (A^-2) = -4=. There exist three branch cuts on the x-plane [-2i,-( v / 3 + 

e)i],[-V3i,V3i], and [(V3 + e)i,2i] before taking the limit xi — > x . These three-branch cuts are joined at the 
locations of x — ±\/3i after taking the limit and they reduce to a single-branch cut [— 2i,2i]. Therefore, we 

can explicitly calculate the values (27V , Ni) as follows: 2N = ^ § Ao T{x)dx = ^ J-Jzi v^+4 ^ x = ^' anc ^ 
N ^^l § Al T ( x ) dx = in I%i vi+4 dx = 1 as we expected. 
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also from the pure U Sp(A) gauge theory (2N — Nf = 4) with the breaking pattern U Sp(A) — > 
USp(2) x 17(1) [15]. 

2) For the second solution USp(Q) — > USp(2) x Z7(2), the matrix model curve y^ 2 = 
t(t + d) 2 — 4?]A 6 has an extra double root at d — — 3eA 2 , e 6 = 1 and it can be factorized as 
yZ 2 = (t- eA 2 ) 2 (t - 4eA 2 ) with rj = e 3 . If we simply put the values (2A" , Aq, N f ) = (2, 2, 2) 
into our general formula (4.1) for the matrix model curve at the singularities, we obtain the 
following results, 

y 2 m = x 2 (x 2 + 3??A 2 ) 2 + 4r? 3 A 6 = (x 2 + ??A 2 ) 2 (x 2 + 4r?A 2 ) , 

where 77 is 6-th root of unity. If we identify this rj with — e in the solution for factorization 
problem, we see the agreement of two approaches. 22 At the intersection points, the vacua 
survive when the M = 2 theory is perturbed by a quadratic (n = 0) superpotential. 

• USp(<6) with N f = 4 

The factorization problem [20] turned out that the matrix model curve has an overall x 2 
factor and is classified as a degenerated case which will be studied later. 

• USp(Q) with N f = 6 

Again the factorization problem [20] turned out the matrix model curve has an overall x 2 
factor and it is classified as a degenerated case which will be studied later. 

As in the SO(N) case we can obtain the coupling constant by using simple replacement of 
(2A" + 2) in pure gauge theory result with (2N — Nf + 2) we can get the coupling constant 
near the singularity without explicit calculation, 

J_ _ d 2 f p (S k ) 1 " d 2 F p (S k ) ( 2N -N f + 2 \ d 2 f p (S k ) 
2m TtJ dS t dS 3 ij Ni^ 1 dSidSi H N t JdS dSr 

where i,j = 1, 2, • • • , n. In the quartic tree-level superpotential case (n = 1), there is only one 
coupling constant and it is given by, 

1 _ in (2N - N f + 2) 2 1 
2« r " "T6(2iV-iV / + 2) 2 log^)' 

where k' = °> \ — ^y. When the condition (2A^o — Nf + 2) = is satisfied, r becomes zero. As 

x l\ x 2~ x o) 

discussed SO(N) case, the asymptotic freedom for USp(2N ) gauge theory breaks down and 
effectively the situation is the same as n = case. Our result gives a consistency with n = 
case. 

22 By using y^ 2 = (x 2 + 4) (x 2 + l) 2 , and B s (x) = x 2 (x 6 + 6a; 4 + 9a; 2 + 2), the function T(x) is given by 
T(x) = —^=. Therefore, we can explicitly calculate the values (2iVo, iVi) as follows: 2N§ = ^ f_ i ^® i+i dx = 
2, and N 1 = ^ [ 2l -r==dx = 2. 
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5 The curve for degenerated U Sp(2N) gauge theory with 
massless flavors 



Again, we want to get a matrix model curve for a degenerated USp{2N) case from the 
degenerated SO(N) results discussed in previous section. From the matrix model curve for 
SO(N) gauge theory, y 2 md = (x 2 + a 2 ) (x 2 + b 2 ) , by simply replacing (a 2 ,b 2 ) of SO(N) with 
(—a 2 , — b 2 ), we obtain the matrix model curve for USp(2N) gauge theory. In other words, we 
replace the parameter m of a superpotential in SO(N) with — m. Therefore, our matrix model 
curve 23 becomes 



9 / 9 \ 2 „ „ K 2 - F 

V m ,d = ( x -m) +AF, m= - , K 



' _p2\2N-N f +2 



(-F) 



2N -N f +2 

(5.1) 



where we use the notation USp(2N) -> USp(2N ) x C/(iVi) with 2iV = 2A" + 2JVi. 

As in SO(N) case there is a special case in which the condition (2 No — Nf + 2) = is 
satisfied and the matrix model curve is given as 

2 4A 4 

ld^=(^ + ^) V 2Nl = l, for (27V - N f + 2) = 0. (5.2) 

Let us demonstrate these general features (5.1) and (5.2) explicitly by comparing them with 
the results from strong-coupling approach obtained in [20] already. Let us consider U Sp(2N) 
gauge theories with Nj(< 2N + 1) flavors where N = 2 or 3. 

• USp(A) with Nf = 2 

In [20], the first kind solution belongs to the degenerated case. To see the equivalence we 
put the values (2A" , A^i, N f ) = (0, 2, 2) into our formula (5.1) (Note 2N = 2A" + 2A^i). However 
in this breaking pattern, since we have the relation (2N — Nf + 2) = 0, the matrix model curve 
(5.2) is given by 

The factorization problem [20] resulted in y^ 2 = (x 2 — « 2 ) 2 =f4A 4 . The gauge group breaks into 
USp(A) — > USp(0) x U{2). One can see the equivalence between two approaches by identifying 
a 2 = -D. 



23 



The relation between W^(x) and ^(a;) can be written as 

V 2 m , d = F 4 (x)=(^-) 2 +4F, N f <2N+l, 



x 



yl 4 = F A {x)=(p^ +4F-4r 2 A 2 , N f = 2N+l. 

In particular, when Nf = 2N+1, the matrix model curve d has an extra contribution from Ax 2 A 2 . Therefore 
in this particular case, the curve should be d = (x 2 — m — 2A 2 ) 2 + 4F and m + 2 A 2 = K ^ F ■ 
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• USp(A) with N f = 3 

The first kind of solution in [20] is the one corresponding to the degenerated case with the 
breaking pattern y^ 2 — (t + D) 2 + AF, D = — yr. The —a 2 corresponds to D here. To see 
the agreement, we put the values (2N , Ni, Nf) = (2, 1, 3) into (5.1) (the breaking pattern with 
(N , Ni, Nj) = (0,2,3) where USp(A) -»• USp(0) x U{2) reproduces the same result) and get 
the matrix model curve, 

2 / 2 \ 2 A 6 F 2 

y m ,d=( x -m) +4F, m = ~F~AS' 

Therefore, there is an agreement. 

• USp(A) with Nf = A 

1) The first kind of solution in [20] is represented as y m = (t + d) + Ar]A 2 t, rf = 1 and 
the breaking pattern is USp(0) x U(2). The a 2 corresponds to — d ± 2A 2 . If we change 
the parametrization for matrix model curve by d — ^ — r/A 2 , it is rewritten as y^ 2 = 

(t + r]A 2 - -^f + AF. To see the agreement, we put the values (2N ,N u Nf) = (0,2, 4) 
into our general formula (5.1) and get the matrix model curve as 

o / 2 \ 2 TjF 2 

y m ,d=[ x -m) +AF, m= — -r]A 

where rf = 1. We see an agreement. 

2) For the other breaking pattern USp(A) — > USp(2) x U(l), we check the equivalence 
by putting the values (2N ,Ni,Nf) = (2,1,4) into (5.1). In this case since the condition 
(2N — Nf + 2) = is satisfied, the matrix model curve (5.2) is given by 

y 2 m4 =(x 2 + D) 2 -AA\ 

Since the results from the factorization problem is ym 2 = (t — a 2 ) 2 — 4A 4 , we see a perfect 
agreement. 

• USp(A) with N f = 5 

1) For the first kind of solution, the matrix model curve can be written as y~m 

AF, D = e (Jpj 3 - eF (|§) 3 where e 3 = I. The D corresponds to a/2 in [20]. We use different 
parametrization and the breaking pattern is USp(A) — > USp(0) x U(2). Putting the values 
(2N , Ni, Nf) = (0, 2, 5) into (5.1) with the footnote 23 we get the matrix model curve as 

/ n2 / F 2 \ * F { A 2 \ 5 

y 2 m4 = (x 2 - m - 2A 2 ) +AF, m = -r, J + - - 2A 2 

with T] 3 = 1. An agreement between two approaches is evident. 

2) By putting the values (2iV , iVi, iV/) = (2,1,5) into (5.1) with the footnote 23 where 
USp(A) -> USp{2) x 17(1), we get 

y 2 m4 =(x 2 -m- 2A 2 ) 2 + AF, m = ^~ 3A '' 
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From the results in [20], the second kind of solution can be represented as y m 2 — (t + D) 2 + 
4A 2 (3A 2 + to). If we introduce F = A 2 (3A 2 + to), we get to = — 3A 2 + ^- and see an agreement. 

• USp(Q) with N f = 2 

The matrix model curve is iff 2 — (t — f ) — 4eA 4 , where e 3 = 1 and the breaking pattern 
is USp(Q) -> USp(0) x C/(3). In this case since the condition (2N - N f + 2) = is satisfied, 
the matrix model curve (5.2) from glueball approach is given as 

n / n N2 4A 4 fi 

Therefore, by identifying (^r,e) with (D,r] 4 ), we can see the agreement between two results. 

• C/Sp(6) with JV/ = 4 

1) For the first kind of the solution from the factorization problem it was given in [20] 
as yff 2 = (t + D) 2 + 4F, D = r\ F ^f 1 ff A , where vf = 1. We use different parametrization to 
see the equivalence easily. By putting the values (2Ao, N±, Nf) = (4,1,4) where USp(6) — > 
USp(A) x U(l) (We also get the same results if we put (2N ,Ni,N f ) = (0,3,4). Therefore, 
we cannot tell which one corresponds to (^ § A T(x)dx, ^ § Ai T{x)dx) on this branch without 
doing these computations), we obtain the matrix model curve as 

2 / 2 \ 2 lT , V^ 4 -xF^/F 

y m ,d = \ x - m ) + 4F ' m = ^F~ r] H^ r 

where rf = 1. Therefore, the results obtained from two approaches are the same. 

2) For the second solution, the matrix model curve is y^ 2 = (t + d) 2 — 4eA 4 , where e 2 = 1 
and the breaking pattern is given by U Sp(6) — > U Sp(2) xU(2). In this case since the condition 
(2 N — Nf + 2) = is satisfied, the matrix model curve (5.2) from glueball approach is given as 

2 / 2 \ 2 4A 4 a 

Therefore, by identifying (d,e) with (D,r) 2 ), we can see the agreement between two results. 
This solution is exactly the same as the degenerated solution with USp{^) — > USp(0) x U{2) 
with Nf = 2: We can expect to have this feature from the addition map. 

• USp(<6) with Nf — 6 

1) The matrix model curve is ym 2 — (t + d) 2 — 4eA 2 t, where e 2 = 1 and the breaking pattern 
is given by USp(6) -> USp{2) x U(2). In this case by putting the values (2A^ , N u Nf) = (2, 2, 6) 
into (5.1), the matrix model curve from glueball approach is given by 

y 2 m, d = ( x2 - m f + 4F > m = %- 7 i^ 2 

where i] is 2-nd root of unity. Therefore, by identifying F = deA 2 — A 4 and r\ = e respectively 
we can see the agreement between two results. These solutions are exactly the same as the 
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degenerated solution with USp(A) — > USp(0) x U(2) we have discussed (Nf = 4): We can 
expect to have this feature from the addition map. 

2) For the breaking pattern USp(Q) -> C/Sp(4) x 17(1), the condition (2JV - N f + 2) = is 
satisfied and the matrix model curve is given by 

which agrees with the one in [20]. The D here corresponds to —si there. We can see this 
solution from the breaking pattern USp(4) — > USp(2) x U(l) through the addition map. 24 

In summary, for the USp(2N) gauge theories with Nf(< 2N + 1) flavors where N = 2,3, 
the solutions (5.1) and (5.2) coincide with the matrix model curve from the strong-coupling 
approach except that we could not check the USp(6) gauge theory with Nf = 6 where there 
exists a breaking pattern USp(6) — > USp(0) x U(3) (in this case, N = 0). Since the SW 
curve for U Sp(2N) — > U Sp(0) x U (N) has less power of x 2 compared with the SW curve for 
USp(2N) — > USp(2N ) x U(Ni),No ^ 0, the factorization problem will be more complicated 
due to the existence of many parameters. Therefore, one expects that our matrix model curve 
(5.1) and (5.2) will predict the solutions for the degenerated USp(2N) gauge theories with 
N f (< 2N + 1) massless flavors where iV > 4 and USp(2N) -> USp(2N ) x U(N 1 ), N ^ 0. 

6 The curve for degenerated SO(N) gauge theory with 
massive flavors 

For the massive case, there exists one extra condition W^(±m/) = where mf(^ 0) is the 
mass of flavor which is different from the m introduced in [15] or section 2 as a parameter of a 
matrix model curve. From this condition, the mass of flavor can be written as 

m ) = ~\ ( x o + x i + x ^) ' W'^x) = x (x 2 - m)) . 

I — 2" ^ S~ 

Equivalently, rrif = ±iy x o+ x i+ x 2 j s p Ure imaginary where the matrix model curve is described 
by Um — Hi=o (^ 2 + x 1 ) with xq < x\ < x 2 . The fact that rrif is pure imaginary implies 
the flavor-dependent part Flavor can be written in terms of the dual period Hi (not n ) and 
other term. Therefore, we encounter a complete different story when the masses of flavors are 
nonzero. The absolute value \mf \ is greater than xq\ xq < |to/|. When we take the singular limit 
x\ — > Xq, the three-branch cuts [—ix 2 , —ixi], [— ix , ix ] and [1x1,1x2] are joined at x — ±ix . 

24 For the last breaking pattern USp(6) — > USp(0) x ?7(3), since the factorization problem was not solved 
exactly in [20], we cannot check the results explicitly. We simply put the results from glucball approach, 

V 2 m , d = (* 2 + + 4F, and m = # _ ElArf = 1. 
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As in massless cases, we want to study two topics, 1) the behavior of n = and n = 1 
singularities and 2) the curve for degenerated case. Although for the degenerated case, we 
obtain a general matrix model curve (in the sense that the curve holds for the generic value of 
N), as we will see below, there exist some difficulties for n = and n — 1 singularities which 
will be discussed in next section. Thus, in this section we concentrate on the degenerated case 
in which the matrix model curve can be written as y 2 m = x 2 (x 2 + a 2 ) (x 2 + b 2 ). According to 
the notation above, we simply put the values (xo, x±, x 2 ) as (0, a,b). There are two-branch 
cuts [—ib, —ia] and [ia, ib\. The mass of flavors can be written as mj = ±i\J as before. 
Therefore, the position of mass is always on the branch cuts: a < |m/| < b. Recall that the 
property of the integral 25 , 

Mo Mo Mo 

/ y m dx = / y m dx = / y m dx, ix < P < ix 2 . 

Jixo JP J 1X2 

In addition, we have to pay attention to the form of J- 'flavor- Naively the contributions from 
3~ flavor is given by (2.1). However, to get the phase factor correctly it should be as follows and 
by simple manipulations it is written as a dual period ITi plus other term: 

Mo f~ m f Mo 

^flavor = N f / y m dx + N f y m dx = 2iriN f U 1 + 2N f / y m dx. 

Jrrif J—Aq JiAg 



Note that we used the fact that y m dx = irilli. In previous paper [15], there was a typo in 
a dual period rix . Therefore, the effective superpotential with massive flavors in degenerated 
branch can be given as 

W cS = 2m [(N - 2)II + (2N, - N f ) U 1 ]—2(N — 2 — N f ) S log f A 



A , 

-2N f I ° y m dx. (6.3) 

JiA 

It is ready to calculate the equations of motion for the fields f and / 2 , in general. When 
we compute the equation of motion for a field F for degenerated case, it is noteworthy to 
observe the following relation under the large A limit where y 2 ^ = x 2 (x 2 + a 2 ) (x 2 + b 2 ^ 



x 2 



2 

'x 2 -m 2 f ) +4F 



2Nf [ A ° i^dx ~ -4N f \ogi. 

JiAct Or 



liAo dF 



25 One can consider a noncompact cycle B\ which has an intersection number one with A\ compact cycle 
and meet at x = rrif = P. Then since the closed loop B\ — B\ does not contain any pole in the y m (x), 
there exists a relation § _~ y rn dx = 0. Therefore, the dual period lb = 7^7 § B y m dx can be written as 

fiAo 1 , 1 rixi 1 \ j t ■ ^ mi 1 1 r*Ao 



fix" y mdx + 2¥I Jill {-Vm)dx which will be X f ix ° y m dx. We perform similar computation for the cycle B\ 
and obtain X jp Aa y m c 
fp° Vmdx = y m dx. 



and obtain X j p Aa y m dx. By changing the upper bound iA into A with an extra term, we get the relation 
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By using previous results given in [15] 26 and this log i-term, we obtain one equation of motion, 

4A 2 4A 2 



(7Vo-2)log 



(a + by 



+ (27V! - N f ) log 



(b 2 



Nf log i = 0. 



(6.4) 



Note the last term for the phase factor which comes from the last term in the formula (6.3). 
Strictly speaking, the coefficient in the last term (6.4) we got from the glueball approach is 2, 
not — 1. It is not clear how this arises but on the other hand, the strong-coupling approach 
implies that the relation (6.4) should be correct. Therefore, after we solve this equation, the 
matrix model curve 27 can be represented as 



4 = (x 2 - m}) 2 + 4F, 



K 2 



m f = 



K 



K = 



(eA 2 ) N - 2 (etA 
( V iVF) 2N i 



2\-N f 



-N f 



N -2 



(6.5) 



where e 2 



1. Recall that m) = - (a 2 + b 2 ) /2 and AF = - (b 2 - a 2 ) 2 /A. The quartic (in 



x) matrix model curve in the breaking pattern SO(N) — > SO(N ) x U(Ni) with Nf massive 
flavors (N — N + 2Ni) depends on N , N±, Nf and a field F. Contrary to the non-degenerated 
case, this formula is always valid (On the non-degenerated case, there exists another condition 
c < where c is given by (7.10) or (7.11) which will be discussed in next section). Note that 
we use a different notation to see the condition W^(±rrif) = easily. When the condition 
(N — 2) = is satisfied, we have to use different formula, as we did before, 



v/ 2 

i>m,d 



(x 2 - 



m 2 ) 



4A 4 
rf 



rj 



2N-i-Nf 



Here we used the equation (6.4) when (Aq — 2) 



for (N - 2) = 0. 

A A 2 

and identified r\ with ± fc2 _ a2 . 



(6.6) 



Note 



mj = -i(a 2 + b 2 ) 



The discussions above are valid to the r = branch studied in [19]. The immediate question 
is how do we represent the matrix model curve for r ^ vacua? As already discussed in [6], 
the nonzero r- vacua can be realized by changing a mass parameter m/. If the singularity 
passes through the first cut (enclosed by the A\ contour), the A\ contour is deformed and 
has transformed to A 1 + C\ — C\. Note that in our present case, the contours A\ and C\[C\) 
correspond to a contour around the branch cut [1x1,1x2} and a contour around m/ on the 
upper(lower) sheet respectively. On the pseudo-confining vacua the singularity is located on 
the second sheet, so we have §g T(z)dz = and ^ §c 1 T(z)dz = 1. After transition we 

4A" 



26 There are two relations for the derivatives of dual periods with respect to a field F: 2m^^- = 2 log | ^ b yi \ , 
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has an extra contribution from 4A 2 (x 2 — m 2 ) . 

- 2 , 

K 2 -F 



and 2?rifii 



27 In particular, when Nf = N — 3, the matrix model curve y 2 ^ d 



Therefore, in this particular case, the curve should be y 2 



7 



2A 2 



- AF and m 2 ■ 



-2A 2 



K 
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have N x — 1 = ^- § Ai T(z)dz. Thus, the effective value of N x is reduced by 1 in this transition. 
Therefore, if we pass through r-singularities, after transition the effective value of N x has been 
reduced by (Ni — r). In addition, the effective number of flavor is also reduced by (Nf — r). 
The contribution from the singularity on the first sheet has different sign from the one on the 
second sheet because the sign of y m is different between the two sheets. Therefore, by adding 
the quantity r from the second sheet, J- 'flavor on the r-branch can be summed and represented 
as follows: 

Mo r~ m f Mo r~ m f 

^flavor = (N f -r) y m dx + (N f - r) y m dx + r (-y m ) dx + r (-y m ) dx 

Jmj J— Ao Jirif J—A-o 

= 2Nf [ ° y m dx = 27niV/ITi + 2N} [ ° y m dx 

J ix2 J iAo 

where we introduce the following notation, 

N = N - 2r, = N f - 2r, N ± = N x - r. (6.7) 

Therefore, by combining the modified N± and the contribution from the flavors, after transition 
the effective superpotential can be given as 



W eff = 2vri 



(N - 2)n + (27V! - N f ) IT] - 2(N - 2 - N f )S log 



-2N f 



Mo 

/ y m dx. 

JiA 



In the log term, we used the fact iV — 2 — Nf = N — 2 — Nf. This effective superpotential for 
nonzero r-vacua takes the same form for r = vacua (6.3) with modified quantities (6.7) and 
therefore includes (6.3). 

Now it is straightforward to obtain the matrix model curve for r ^ vacua. The general 
matrix model curve 28 can be given by starting with (6.5) and changing the values (N, N±, Nf) 
appearing in (6.5) into (N, N ly Nf) where they are defined as in (6.7) : 



y 2 m>d = (x 2 -mj) 2 + 4F, m J = K 



(eA 2 ) jV - 2 (6^A 2 )- jV ' ' 



JV -2 

(6.8) 



where e 2 = rf = 1 and moreover when the iV is equal to 2, as we did before, the matrix model 
curve with modified quantities (6.7) becomes 

yl 4 =(x 2 -m}) 2 -^ Tr>-»'=i»', for (JV -2) = 0. (6.9) 



28 When Nf = N - 3, the matrix model curve should be y 2 m d = (z 2 - m 2 - 2A 2 ) + 4F and mj + 2A 2 = 



k ■ 
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Therefore both equations (6.8) and (6.9) are the most general expressions including (6.5) and 
(6.6). To demonstrate these general solutions, let us consider some explicit examples, SO(N) 
theories with Nf(< N — 3) flavors for TV = 4, 5, or 6. For SO (6) gauge theories, we only consider 
even number of flavors as we did in [19]. 

• 5*0(4) with Nf — 1: Non-baryonic r = branch 
1) 50(4) -> 50(2) x [7(1) 

The matrix model curve derived by the factorization problem of degenerated case was given 
in [19]. After the factorization (Note that the F 6 (x) has an overall factor x 2 in the non- 
degenerated case) we have obtained F 4 (x) = (x 2 — A) 2 — 4A 2 (x 2 — m 2 ^ and from the relation 
between W%(x) and F 4 (x) one can read off the W^(x). By the condition W^(±mf) = 0, we 
see A = m 2 and the matrix model curve from strong-coupling approach can be written as 
yZ 2 = (x 2 - m} - 2A 2 f -4A 4 . 

Now let us see this observation from the glueball approach. If we substitute the values 
(iVo, Ni, Nf) = (2, 1, 1) where N — N + 2iV\ into our general formula (6.6) (In this case the 
condition (N — 2) = is satisfied and also see the footnote 27), we obtain the matrix model 
curve on the degenerated case, 

y 2 m4 =(x 2 -m 2 -2A 2 ) 2 + 4A* 

with the phase factor rj — i. Therefore, we do not see the agreement:^ 2 ^ Umd- Let 
us go back a derivation for (6.6). According to the strong-coupling approach, the curve is 
characterized by y^ 2 = (x 2 — m 2 ^ (x 2 — m 2 — 4A 2 ^J . Then a 2 in the matrix model curve 
Vrnd = (x 2 + ° 2 ) ( x2 + corresponds to — m 2 — 4A 2 and b 2 corresponds to —m 2 . So in this 
particular case, the previous condition a < |m/| < b does not hold. In other words, the glueball 
approach is not allowed. 29 

• 50(5) with Nf = 1: Non-baryonic r = branch 

To see the equivalence between two approaches easily, we solve the factorization again by 
using the following parametrization, Pl(x) — 4a; 2 A 4 (x 2 — m 2 ^ = H%(x) [(x 2 — b) 2 + 4F]. Note 
AF = c in the notation of [19]. From the relationship between F 4 (x) and W^x), we obtain 
b = m 2 and additionally there exists one equation corresponding to (5.15) of [19]. However 
for the present purpose, we reexpress it in terms of F by writing a as a function of F or c, 

29 The first kind of solution in [19] where SO{A) — > U(2) has the matrix model curve as y^' '■ = (x 2 -D) +AF 
where D = —6/2 and AF = c — b 2 /A in the notation of [19]. In order to get these results, we used two 
relations: One for c and the other is a 2 (a + 9A 2 ) + 3a (a + 4A 2 ) m 2 + (3a + 4A 2 ) m 4 - + m f — together with 
b = — 2m 2 — 4A 2 . By using these two relations, one can get an equation for m 2 or D in terms of F or c. To see the 
equivalence between two approaches (glueball approach and strong-coupling approach), we put (No, Ni,Nf) = 

(0, 2, 1) into the general formula (6.5) and get the matrix model curve d = (x 2 — m 2 — 2A 2 ^j + AF, and 

3 

m 2 = —^^ — h F—^-rf + 2 A 2 where r/ is 4-th root of unity. This leads to a disagreement. 
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F 4 -2F 2 A 8 + A 16 -FA 8 m 4 = 0. Solving this equation, we can get mj as follows: mj = ±^=£r- 
In this case, the matrix model curve is written as y^ 2 = (x 2 — mj^j + AF. 

On the other hand, if we put the values (N , N u N f ) = (3, 1, 1) where 50(5) -> SO (3) x[7(l) 
into (6.5) we obtain the matrix model curve, 

9 / 9 ?\ 2 9 (A 8 -F 2 ) 

^ = (* 2 -™?) + 4i? > m / = ± aV F • 

Therefore, there is a perfect agreement. Although there is a breaking pattern 5*0(5) — > 50(1) x 
C/(2), the corresponding matrix model curve does not exist. 

• 50(5) with iVj = 2: Non-baryonic r = 1 branch 

From the solution for the factorization, we have obtained y^ 2 = (x 2 — s) 2 — Ax 2 A 2 where we 
wrote the characteristic function Pi{x) in [19] as (x 2 — s). Taking into account the condition 
W^(±mf) = 0, we get mj = s. If we introduce a new notation AF = — 4A 4 — As A 2 , the mass of 
flavor can be rewritten as mj = —A 2 — and the curve becomes y^ 2 = (x 2 — m 2 — 2A 2 ) + 
AF. Since the breaking pattern of this solution is 50(5) — > 50(1) x £7(2), by putting 
(TV, Ni, Nf, N ) = (3, 1, 0, 1) into our general formula (6.8) (we took e = —1) with the footnote 
28, we obtain the matrix model curve, 

y 2 m4 = (x 2 -m 2 - 2A 2 ) 2 + AF, m) = -A 2 - ^. 

Therefore, we see the agreement between two approaches. 

• 50(5) with Nf = 2: Non-baryonic r = branch 

By using the same parametrization as [19] with F 4 (x) = (x 4 + Abx 2 + 2c) we can represent 
c as a function of a, b and rrif and there is an extra equation characterized by, a 2 + mj + 
2a (2A 2 + m 2 f) = 0. To see the equivalence between two approaches let us introduce a new 
notation AF = c—Ab 2 and take into account the condition W^(±m/) = 0. Finally, by taking 2b = 
—m 2 — 2A 2 and solving the two equations (the relation for c and the above equation) we get the 
mass m 2 = — 3A 2 +^-. The matrix model curve can be written as y~„ 2 = (x 2 — m 2 — 2A 2 ) +AF. 
If we put (N ,N u N f ) = (1,2,2) where there is 50(5) -> 50(1) x [7(2) into (6.5) (we took 
e = — 1) with the footnote 27, we obtain the matrix model curve 

y 2 nA = (x 2 -m 2 f - 2A 2 ) 2 + AF, m 2 = ^ - 3A 2 . 

We see an agreement. In this case, for (mj + 2 A 2 ) + AF = 0, there is a Chebyshev branch 
with 50(5) -> 50(5). 

• 50(6) with Nf = 1: Non-baryonic r = branch 

If we rewrite c in the results of [19] as c = AF the matrix model curve for degenerated case 
can be represented as y^ 2 = (x 2 — mfj + AF with one equation, F 6 — 2F 3 A 12 — AF 2 mjA 12 — 
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Fm 8 fA 12 + A 24 = which is the last equation of (5.18) in [19]. Solving this constraint, we can 



obtain m 2 f as follows: m 2 f = ±- or ±- r 

To see the agreement with our general formula we simply put (No, JVi, iV/) = (4, 1, 1) where 
there exists SO (6) — >■ 5*0(4) x £/(l) into (6.5) we can get the matrix model curve 

2 / 2 2\ 2 2 (A 6 -^) (A 6 + F§) 

2/m d = ( x ~ m f) + 4 ^ m f = ± A 1 ±i- 1 -■ 

This implies a perfect agreement. For m 4 . + 4F = 0, a Chebyshev branch 50(6) — > 50(6) 
appears. 

• 50(6) with Nf = 2: Non-baryonic r = branch 

There are two types of solutions in [19], one has a t factor in ^(t) which can be written as 
tH\(i) and the other has no t factor in H2(t). 

1) 50(6) -> 50(2) x £/(2) 

First, let us consider the former case. The matrix model curve given in [19] is y^ 2 = 
(t - m 2 ff + 4A 4 by inserting a = mj + 2r?A 2 into F 4 (x). By putting (jV , A r i, N f ) = (2, 2, 2) 
into our general formula (6.6) we obtain the matrix model curve from the viewpoint of glueball 
approach, 

^=(- 2 "^) 2 + 4A 4 . 

This leads to an agreement between two approaches. When Z} 2 + 4A 4 = 0, there is a Chebyshev 
branch 50(6) -> 50(6). 

2) 50(6) -> f/(3) 

Secondly, we move to the latter case. The solutions were displayed in equation (5.22) of 
[19]. If we use a new notation F = — (b + A 2 ) 2 , the matrix model curve can be rewritten as 
y^ 2 = (x 2 — m 2 fj +AF with m 2 = ±i ^—jt - . To see the equivalence with our general formula, 
we put the breaking pattern on this solution, namely 50(6) — > U(3) into the formula (6.8) and 
get the matrix model curve 

yl, d = {- 2 -m}) 2 + AF, m ) = ±^d^±. 

One can see an agreement between two approaches. 

• 50(6) with Nf = 2: Non-baryonic r = 1 branch 
1) 50(6) -> 17(3) 

The first kind of solution studied in [19] gives a relation for mass m 2 = a — 2r]A 2 . There was 
a typo in the equation of [19]. If we introduce a new notation AF = Aar/A 2 — 4A 4 , the mass can 
be rewritten as m 2 = r] ^ — A 2 ^. Since this solution has the breaking pattern 50(6) — > U(3) 
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by putting (N, N ± , Nf, N ) = (4,2,0,0) into our general formula (6.8) we obtain the matrix 
model curve 

V 2 m,d = (x 2 - m}) 2 + AF, m) = -e - A 2 ) 

where e is 2-th root of unity. This leads to an agreement. 
2) SO{6) -> 50(2) x [7(2) 

The other solution has the matrix model curve y^ 2 = {x 2 — m 2 ) — 4A 4 . Since the condition 
(N — 2) = is satisfied, the matrix model curve from (6.9) can be rewritten as 

A perfect agreement occurs. 

• 5*0(6) with Nf = 3: Non-baryonic r = branch 

The solution for the breaking pattern SO (2) x U(2) can be represented as follows: y^ 2 = 
(x 2 — m 2 — 2 A 2 ) — 4A 4 . We use F = A 2 [p + m 2 + A 2 ) . From the relationship between F 4 (x) 
and W%(x), we have D = —m 2 — 2A 2 implying F = —A 4 . Putting this relation we get the above 
matrix model curve. 

On the other hand, by substituting the values (No, Ni, Nf) = (2,2,3) where 50(6) — > 
50(2) x U(2) into (6.6) we obtain the matrix model curve as 

*4,= (* 2 + £) 2 + 4A 4 . 

However, as we have discussed in 50(4) with Nf = 1 case, the glueball approach given by (6.6) 
does not include this case also: The location of the flavors is different from the region belonging 
to the branch cut. For D 2 - 4A 4 = 0, 50(6) -> 50(6) Chebyshev branch occurs. 

• 50(6) with Nf = 3: Non-baryonic r = 1 branch 

In [19], it was discussed by using the addition map that this branch was the same as the 
one for 50(4) with Nf = 1 case. Since our general formula also has the same structure, the 
glueball approach here reproduces the same results from the factorization problem there. That 
is, for the breaking pattern, 50(6) — > U(3), the values (N , N ± , Nf) = (0,2,1) correspond 
to (N , Ni, Nf) = (0,2,1) for 50(4) theory with N f = 1. Moreover, on the breaking pattern 
50(6) -> 50(2) xU(2), the values (N , Ni,N~f) = (2, 1, 1) correspond to (N ,N u N f ) = (2,1,1) 
of 50(4) theory with N f = 1. 

In summary, for the SO(N) gauge theories with Nf(< N — 3) flavors where N = 4, 5 or 6, 
we have checked that the solutions (6.5), (6.6), (6.8), and (6.9) coincide with the matrix model 
curve from the strong-coupling approach except that we could not check the 50(4) gauge theory 
with Nf = 1 where there exists a breaking pattern 50(4) — > U(2) (in this case, = 0) and 
50(6) with Nf = 3 where there exists a breaking pattern 50(6) — > U(3) (N = 0). Therefore, 
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at least one expects that our matrix model curve (6.5), (6.6), (6.8), and (6.9) will predict the 
solutions for the degenerated SO(N) gauge theories with Nf(< N — 3) massive flavors where 
N >7 and SO(N) -> SO(N ) x JV ^ 0. 30 



7 Discussion on the n = and n = 1 singularity with 
massive flavors 

To obtain a general matrix model curve for n = and n = 1 singularities for our gauge 
theories, we have to describe the behavior of the matrix model curve y 2 m = Ylj =0 (x 2 + x 2 ) in 
the limit x\ — > x u . Contrary to the degenerated case we have studied in previous section, there 
exist two possibilities for the position of the mass of flavors. One is x < \rrif\ < x 2 and the 
other is x 2 < < \/f x 2- In the former case, since the position of flavor mass is on the 
branch cut [1x1,1x2], J 7 flavor can be written in terms of a dual period ITi. Then by following 
the derivation in previous section 31 we obtain the matrix model curve near the singularity on 
the non-degenerated case as follows: 

y 2 m = x 2 (x 2 -m 2 f ) 2 + (f 2 )x 2 + (f )=[x 2 ±2r ] A 2 (l + c)] 2 (x 2 ±4r ] A 2 ), 

where the fields, parameter of a superpotential, and glueball field are given by 



m 



(f 2 ) = VA 4 (l + 2c), 
J = ±2 V A 2 (2 + c), 

n(2N 1 - Nf) 
N-N}-2 



(/ ) = ±16^A b (l + c) z , 
(S) = - V 2 A 4 (1 + 2c) , 



(7.10) 



c = cos 
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N—Nf—2 



i N f 



together with (6.7). We can see that the solutions for the roots x\ = ±4i]A 2 and x^ = 
±2i]A 2 (1 + c) satisfy the relation 2xq < x\, implying that these solutions are consistent with 
the condition we have assumed before (Note |m/| = J (2xq -\-x 2 )/2 < x 2 ). Note that c should 



30 Using this solution from the glueball approach, one obtains the matrix model curve for arbitrary N. For 
example, for SO(7) with Nf = 1 (non-baryonic r = branch), when we consider the breaking pattern 5*0(3) x 

U(2), by putting (No, Ni, Nf) = (3,2, 1) into our general formula (6.5) we predict the matrix model curve as, 

/ \2 (A 16 —F 4 ') 

Vmd = [x 2 ~ m / ) and m 2 f = ± F fp AS ■ F° r SO (7) with Nf = 2 (non-baryonic r = 1 branch), after we 

put (N, Ni,Nf, N ) — (5, 1, 0, 3) into (6.8) we obtain the matrix model curve, V^ n ^ = (x 2 — m^j + 4F, and 

mj = — + (In this case, we took e = —1). Moreover, for SO(7) with Nf = 3 (non-baryonic r = 1 branch), 
the factorization problem becomes the same as the one for SO(5) with Nf = 1 case through the addition map. 

31 Note that when we compute the equation of motion for fo, since the extra piece from the flavor part, an 
integral ^f^dx, behaves like C(Aq 2 ), it does not contribute. Then the equation of motion for fo looks 



similar to the pure case [15] and it is given by 



cos 



ir(2Ni-N f ) 
N-N,-2 



and the equation of motion for fa 



leads to (n - N f - 2) log|^| + Nflogi = 0. After solving this, one gets x\ = ±4r/A 2 , and V N - N f- 2 = i N f . 
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be negative or zero but c ^ -1. When c = — 1, then this will be the vacuum in the degenerated 
case. The constraint \m,f\ < x 2 is too restrictive and most of the examples from factorization 
problem [19, 20] are beyond this criterion. 

There exists one more possibility, x 2 < \m,f\ < \f\x 2 . In this case, since Ffiavor cannot 
be written in terms of a dual period III, we cannot obtain a general matrix model curve. 
Namely, the integral Tfi avar can be represented by an elementary function but when we solve 
the equations to obtain the roots xq and x 2 , a general solution for the matrix model curve is 
not apparent. 

In summary, at n = and n = 1 singularity we obtain some partial results for matrix 
model curve (the curve does not exist for all the range of N, N 1: and Nf but for some particular 
values) in which the condition, c is negative or zero (c ^ —1), should be satisfied. To see the 
correctness of the result let us study an example satisfying this condition only. 32 

• 5*0(5) with Nf = 1: Non-baryonic r = branch 

The matrix model curve derived by the factorization was given in [19] as (There is a typo 
in the presentation of [19]) y^ 2 = x 2 (x 2 — m 2 ) — 4A 4 (x 2 — rnfj. This matrix model curve 
has an extra double root at mj = ArjiA 2 where rf = 1 and have the following factorized 
form, {jm 2 = {x 2 ± 2ih 2 ) 2 {x 2 ± 4iA 2 ). To see the equivalence of two approaches we substitute 
(N ,N u N f ) = (3,1,1) where SO(5) -> SO(3) x 17(1) and N = N + 2N t into our general 
formula (7.10) and then obtain the matrix model curve: 

y 2 m = (x 2 ± 2r l A 2 f (x 2 ± 4r?A 2 ) 

where rf — %. Therefore for pure imaginary rj = —i, we see the agreement. 

For USp{2N) already discussed in section 4, we have to pay attention to the signs 

in the matrix model curve. To get the right results for USp(2N) theory from the results for 
SO(N) theory, we simply replace x^ — * —x\,x\ — > — x 2 ,, and (JV — 2) — > (2A^ + 2) respectively. 
With this in mind, the matrix model can be written as 

Vl = x 2 (x 2 -m 2 f ) 2 + (f 2 )x 2 + (f )=[x 2 ±2r ] A 2 (l + c)] 2 (x 2 ±4r ] A 2 ), 

where the fields, parameter of a superpotential, and glueball field are given by 

tfr> = 4t? 2 A 4 (1 + 2c) , (/o) = ±16^ 3 A 6 (1 + cf , 

32 For 5*0(5) with Nf = 2 (non-baryonic r — branch), the matrix model curve was given in [19] as 
jMi 2 = (x 2 - m 2 f ) (a; 2 -4A 2 ). To see whether there exists the equivalence of two approaches, we sub- 
stitute (No, N\, Nf) = (3,1,2) into our general formula (7.10) and then obtain the matrix model curve: 
y 2 n = (x 2 ± 2rjA 2 ) 2 (x 2 ± ArjA 2 ) where rj = —1 with m 2 = ±4r]A 2 . However, in this case, c = 1. There- 
fore, we cannot see an agreement between two approaches. That is, the matrix model curve exists only from 
the strong-coupling approach. One can easily see that for SO(6) with Nf = 1 (non-baryonic r = branch), 
the c value is 1/2, for SO (6) with Nf = 2 (non-baryonic r = branch), c = 1, and for SO (6) with Nf = 3 
(non-baryonic r = branch) the c value is —1. Most of the examples studied in [19] belong to the case where c 
is positive or —1. 
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m) = ±2r]A 2 (2 + c), (S) = -rfA A (1 + 2c) , (7.11) 



cos 



7r(2iVi - Nf) 
2iV-iV/ + 2. 



2N-N f +2 = 



with the modified quantities are given by 2N = 2N — 2r, Nf = Nf — 2r, and 2Ni = 2Ni — 2r. In 
addition, we can see that the solutions for different parametrization are given by x\ = ±4t]A 2 , 
and Xq = ±2i]A 2 (1 + c) which satisfy the relation 2x1 — x \-> implying that these solutions are 
consistent with the condition we have assumed before. Note that c defined in (7.11) should be 
negative or zero (c 7^ —1). 33 Let us consider one example. 
• USp{4) with Nf = 2: Non-baryonic r = 1 branch 

The matrix model curve can be written as = t 3 — 2m 2 ft 2 + 4A 4 + m^j t + 4A 4 m 2 
with a = —2m 2 . Requiring an extra double root it becomes 2 = (t±2iA 2 ) 2 (t±MA 2 ) 
at m 2 = =)=4iA 2 . To see the equivalence we put the values (2N,2Ni,Nf) = (2,2,0) where 
U Sp(4) — > USp(0) x U{2) into (7.11) we obtain the matrix model 

y 2 m = (x 2 ± 2r^A 2 ) 2 (x 2 ± 4r^A 2 ) , rf = l. 

For the cases of pure imaginary r\ = ±i, the matrix model curve from the glueball approach 
coincides with the curve from strong-coupling approach. 
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